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1 Introduction 



We are interested in the probabilistic representation of the solution to a 
porous media type equation given by 



in the sense of distributions, where uq is an initial probability density. We 
look for a solution of (jl.ip with time evolution in 

We always make the general following assumption on [3. 

Assumption 1.1 • /3 :M— >M is monotone increasing. 

• \(3(u)\ < const |u|, u G M. 

• There is A > such that ((3 =p \id)(x) — > =Foo when x — > =Foo where 
id(x) = x. 



Remark 1.2 7. Since (3 is monotone, ( fi. i|) implies that (3(u) = $ 2 (tt)u, 
$ feeing a non-negative bounded Borel function. 

2. (3(0) = and /? is continuous at zero. 

We recall that when = |M|ii m_1 , m > 1, (jl.ip is nothing else but the 
classical porous media equation. 

One of our final targets is to consider $ as continuous except for a possible 
jump at one positive point, say e c > 0. A typical example is 



H being the Heaviside function. 

The analysis of (|l-lj) and its probabilistic representation can be done in the 
framework of monotone partial differential equations (PDE) allowing multi- 
valued functions and will be discussed in detail in the main body of the 
paper. This extension is necessary, among other things, to allow the graph 




\dlMu)), t€ [0,oo [ 

uq(x), 



(1.1) 



• n G (L l r\L°°)(R). 



$(n) = H(u - e c ), 



(1.2) 
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associated with (3 to be a maximal monotone graph. We refer to Assumption 
13.11 below. In this introduction, for simplicity, we restrict our presentation 
to the single-valued case. 

Definition 1.3 We will say that equation U.l\) or (3 is non-degenerate if 

there is a constant Co > such that <3? > cq. 

Of course, if is as in (|1.2I) . then (3 in is not non-degenerate. In order to 
have f3 to be non-degenerate, one needs to add a positive constant to it. 

Several contributions were made in this framework starting from [10] for 
existence, [IB] for uniqueness in the case of bounded solutions and [TT] for 
continuous dependence on the coefficients. The authors consider the case 
where (3 is continuous, even if their arguments allow some extensions for the 
discontinuous case. 

As mentioned in the abstract, the first motivation of this paper was to dis- 
cuss a continuous time model of self-organized criticality (SOC), which are 
described by equations of type (|1.2p . 

SOC is a property of dynamical systems which have a critical point as an 
attractor, see [3] for a significant monograph on the subject. SOC is typ- 
ically observed in slowly-driven out-of-equilibrium systems with threshold 
dynamics relaxing through a hierarchy of avalanches of all sizes. We, in 
particular, refer to the interesting physical paper [1J. The latter makes ref- 
erence to a system whose evolution is similar to the evolution of a "snow 
layer" under the influence of an "avalanche effect" which starts when the 
top of the layer attains a critical value e c . Adding a stochastic noise should 
describe other contingent effects. For instance, an additive perturbation by 
noise could describe the regular effect of "snow falling". In Bantay et al. 
([!]) it was proposed to describe this phenomenon by a singular diffusion 
involving precisely a coefficient of the type (|1.2p . 

In the absence of noise the density u(t,x),t > 0,x £ K of this diffusion is 
formally described by (jl.l|) and (3{u) = <&(-u) 2 u where $ is given by (jl.2p . 

Such a discontinuous monotone (3 has to be considered as a multivalued map 
in order to apply monotonicity methods. 
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The singular non-linear diffusion equation (jl.ip models the macroscopic phe- 
nomenon for which we try to give a microscopic probabilistic representation, 
via a non-linear stochastic differential equation (NLSDE) modeling the evo- 
lution of a single point on the layer. 

Even if the irregular diffusion equation (|l.ip can be shown to be well-posed, 
up to now we can only prove existence (but not yet uniqueness) of solutions to 
the corresponding NLSDE. On the other hand if <3? > cq > 0, then uniqueness 
can be proved. For our applications, this will solve the case <&(u) = H (x — 
e c ) + e for some positive e. The main novelty with respect to the literature 
is the fact that $ can be irregular with jumps. 

To the best of our knowledge the first author who considered a probabilistic 
representation (of the type studied in this paper) for the solutions of a non- 
linear deterministic PDE was McKean [25], particularly in relation with the 
so called propagation of chaos. In his case, however, the coefficients were 
smooth. From then on the literature steadily grew and nowadays there is a 
vast amount of contributions to the subject, especially when the non-linearity 
is in the first order part, as e.g. in Burgers equation. We refer the reader to 
the excellent survey papers [52] and [21] . 

A probabilistic interpretation of when (5{u) = \u\u m ~ l ,m > 1, was 

provided for instance in [9]. For the same /3, though the method could be 
adapted to the case where f3 is Lipschitz, in [22] the author has studied the 
evolution equation (jl.ip when the initial condition and the evolution takes 
values in the class of probability distribution functions on BL Therefore, 
instead of an evolution equation in L 1 (M), he considers a state space of 
functions vanishing at — oo and with value 1 at +oo. He studies both the 
probabilistic representation and propagation of chaos. 

Let us now describe the principle of the mentioned probabilistic representa- 
tion. The stochastic differential equation (in the weak sense), rendering the 
probabilistic representation, is given by the following (random) non-linear 
diffusion: 



where W is a classical Brownian motion. The solution of that equation may 




density (Y t ) 
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be visualised as a continuous process Y on some filtered probability space 
(Q,J r ,(J-'t)t>o,P) equipped with a Brownian motion W. By looking at a 
properly chosen version, we can and shall assume that Y : [0, T] x Q — > IR + 
is B([0,T\) (g> ^"-measurable. Of course, we can only have (weak) uniqueness 
for (|1.3p if we fix the initial distribution, i.e. we have to fix the distribution 
(density) uq of Yq. 

The connection with is then given by the following result. 

Theorem 1.4 (%) Let us assume the existence of a solution Y for hi. 3^ . 
Then «:[0,T]xl-t R + provides a solution in the sense of distribu- 
tions of hl.l\) with uo := u(0, •). 

(ii) Let u be a solution of U.l\) in the sense of distributions and let Y solve 
the first equation in hl.3\) with law density v(t,-) and initial law density 
uq = u(0, •). Then 

d t v = \dl x (&{u)v), (1.4) 

in the sense of distributions. In particular, if v is the unique solution 
of jl-4\), with v (0, •) = uq, then v = u. 



Proof. Let ip £ Cp°(R), Y be a solution to the first line of (II. 3p such that 
v(t, •) is the law density Yj, for positive t. We apply Ito's formula to <f(Y), 
to obtain 

tp(Y t ) = <p(Y ) + I y'{Y s )${u{s, Y s )) dW s + \ f <p" (Y S )<S> 2 (u(s , Y.)) ds 
Jo 1 Jo 

Taking expectation we obtain 



1 



y{y)v{t,y)dy = I (p(y)u (y)dy + | / ds \ tp" (y)<S> 2 (u(s,y))v(s,y) dy. 



o 



Now both assertions (i) and (ii) follow. 



Remark 1.5 An immediate consequence of the probabilistic representation 
of a solution of hl.l\) is its positivity at any time. Also the property that the 
initial condition is of mass 1 is in this case conserved. 
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The main purpose of this paper is to show existence and uniqueness in law 
of the probabilistic representation equation (11 .3f) . in the case that (3 is non- 
degenerate and not necessarily continuous. In addition, we prove existence 
for (jl.3h . in some degenerate cases under certain conditions, see Subsection 



Let us now briefly explain the points that we are able to treat and the diffi- 
culties which naturally appear in regard to the probabilistic representation. 

For simplicity we do this for (5 being single- valued (and) continuous. How- 
ever, with some technical complications this generalizes to the multi-valued 
case, as spelled out in the subsequent sections. 

1. Monotonicity methods allow us to show existence and uniqueness of 
solutions to (jl.ip in the sense of distributions under the assumption 
that (5 is monotone, that there exists A > with {j3 + Xid) (R) = R and 
that (3 is continuous at zero, see Proposition 13.41 below. We emphasize 
that for uniqueness no surjectivity of /3 + Xid is required, see Remark 
13.61 below. 

2. Let a : [0,T] x 1R — > R be a strictly positive bounded Borel function. 
Let .M(M) be the set of all signed measures on M with finite total 
variation. We prove uniqueness of solutions of 



as an evolution problem in M (M) , at least under an additional assump- 
tion (A), see Theorem 13.81 below. 

3. If /3 is non-degenerate, we can construct a unique (weak) solution Y 
to the non-linear SDE corresponding to (jl.3p . for any initial bounded 
probability density uq on M, see Theorem 14. 41 below. For this construc- 
tion, items 1. and 2. above are used in a crucial way. 

4. Suppose (5 possibly degenerate. We fix a bounded probability density 
uq. We set <I> e = $ + e and consider the weak solution Y £ of 



K2[ 



d t v 
v(0, x) 



9xx( av ) 

u (x), 



(1.5) 




(1.6) 
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where u £ (t, •) is the law of Yf,t > 0, and Yq is distributed according 
to uo(x)dx. The sequence of laws of the processes (Y £ ) are tight. 
However, the limiting processes of convergent subsequences may in 
general not solve the SDE 



However, under some additional assumptions, see Properties 14.81 and 
14.101 below, it will be the case. The analysis of the degenerate case 
in greater generality (including case (|1.2p ) will be the subject of the 
forthcoming paper [7]. 

In this paper, we proceed as follows. Section 2 is devoted to preliminaries 
about elliptic PDEs satisfying monotonicity conditions. 

In Section 3, we first state a general existence and uniqueness result (Propo- 
sition [33]) for equation (jl.ip and provide its proof, see item 1. above. The 
rest of Section 3 is devoted to the study of the uniqueness of a determinis- 
tic, time inhomogeneous singular linear equation with evolution in the space 
of probabilities on M. This will be applied for studying the uniqueness of 
equation (jl.5p in item 2. above. This is only possible in the non-degenerate 
case, see Theorem 13. H\ if (3 is not non-degenerate we give a counterexample 
in Remark 13.111 

Section 4 is devoted to the probabilistic representation (jl.3p . In particu- 
lar, in the non-degenerate (however not smooth) case, Theorem 14.41 gives 
existence and uniqueness of the non-linear diffusion (|1.3|) which represents 
probabilistically (jl.ip . In the degenerate case, Proposition 14.121 gives an 
existence result. 

Finally, we would like to mention that, in order to keep this paper self- 
contained and make it accessible to a larger audience, we include the analytic 
background material and necessary (through standard) definitions. Likewise, 
we tried to explain all details on the analytic delicate and quite technical 
parts of the paper which form the backbone of the proofs for our main result. 




(1.7) 
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2 Preliminaries 



We start with some basic analytical framework. 

If / : R — > R is a bounded function we will set ||/||oo := sup^eR 1/(^)1- By 
C&(R) we denote the space of bounded continuous real functions and by S (R) 
the space of rapidly decreasing infinitely differentiable functions 92 : R —* R, 
by S' (R) its dual (the space of tempered distributions). 

Let K £ be the Green function of e — A, that is the kernel of the operator 
(e - A)" 1 : L 2 (R) L 2 (R). So, for all v? 6 L 2 (R), we have 

B e (tp) := (e - A)-V(x) = / lf c (a: - y) tp(y)dy. (2.8) 



The next lemma provides us with an explicit expression of the kernel function 



K e . 



Lemma 2.1 



K e (x) = —=e-^ x \x€R. (2.9) 
zy £ 



Proof. From Def. 6.27 in [30], we get 

1 r° _i _m 2 



K *& = 77^ml t-^e-^r-^dt (2.10) 
(47r) i/2 Jo 

The result follows by standard calculus. ■ 

Clearly, if <p E C 2 (R) P| 5'(R), then (e — A) ip coincides with the classical 
associated PDE operator. 

Lemma 2.2 Let £ > 0, m £ M(M). There is a unique solution v e € 

c 6 (M)n(n P >i^w) of 

£v £ — Av £ = m (2-11) 
in the sense of distributions given by 

v £ (x) := [ K £ (x - y)dm(y), x e R. (2.12) 



Moreover it fulfills 



su PV ^k(x)|<^P^, (2.13) 

X ^ 
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where II m\\ var denotes the total variation norm. In addition, the derivative 
v' £ has a bounded cadlag version which is locally of bounded variation. 

In the sequel, in analogy with (|2.8p . that solution will be denoted by B £ m. 

Proof. Uniqueness follows from an obvious application of Fourier trans- 
form. In fact, it holds even in 5'(R). 

v e given by (|2.12p . clearly satisfies (|2.1ip in the sense of distributions. By 
Lebesgue's dominated convergence theorem and because K £ is a bounded 
continuous function, it follows that v £ £ Cf>(M). 

By Lemma 12.11 we have 

sup|w e (x)| < — ^||m|| var , (2.14) 

x £ 

By Fubini's theorem and (|2.9p . it follows that v £ G L 1 (M). Hence v £ E 
L P (M), Vp > 1, because v £ is bounded. 

Since v £ equals ev £ — m in the sense of distributions, after integration, we 
can show that 

v 'e( x ) = £ \ Ve(y)dy-m(]-oo,x]), 

J — oo 

for dx-a.e. x £ M. In particular, v £ has a bounded cadlag version which is 
locally of bounded variation and 

iKlloo < e||«|Ui(k) + ||m|| var . 



We now recall some basic notions from the analysis of monotone operators. 
More information can also be found for instance in [29]. See also [6| 115]. 

Let £ be a general Banach space. 

One of the most basic notions of this paper is the one of a multivalued 
function (graph). A multivalued function (graph) (3 on E will be a subset 
of E x E. It can be seen, either as a family of couples (e, /), e, / S E and 
we will write / S (3{e) or as a function (3 : E — > ~P(E). 

We start with a definition in the case E = M. 
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Definition 2.3 A multivalued function (3 defined on R with values in subsets 
o/R is said to be monotone if (xi — X2){yi — y.2) > for all x\,X2 G R, 
Vi G (3(xi),i = 1,2. 

VFe say i/iat /? is maximal monotone i/ it is monotone and there exists 
A > such that (5 + Xid is surjective, i. e. 

11(13 + Xid) := (J {p[x) + Xx) = R. 

We recall that one motivation of this paper is the case where j3(u) = H (u — 
e c )u. 

Let us consider a monotone function ip. Then all the discontinuities are of 
jump type. At every discontinuity point x of tjj, it is possible to complete ijj, 
producing a multi-valued function, by setting ip(x) = [tp{x— ), ip(x+)]. 

Since ip is a monotone function, the corresponding multivalued function will 
be, of course, also monotone. 

Now we come back to the case of our general Banach space E with norm 
|| • || . An operator T : E — > E is said to be a contraction if it is Lipschitz 
of norm less or equal to 1 and T(0) = 0. 

Definition 2.4 A map A : E — > E, or more generally a multivalued map 
A : E — > 'P(-E') is said to 6e accretive if for all /i,/2, 51,52 £ ^ sitc/i 
<7i G A/j,i = 1,2, we have 

II/1-/2II < ||/i-/ 2 + A( 5l - 52 )||, 

/or any A > 0. 

This is equivalent to saying the following: for any A > 0, (1 + XA)^ 1 is a 
contraction for any A > on Rg(I + XA). We remark that a contraction is 
necessarily single-valued. 

Remark 2.5 Suppose that E is a Hilbert space equipped with the scalar 
product ( ,)h- Then A is accretive if and only if it is monotone i.e. 
(/1 - h, 9i ~ 92)h > for any f u f 2 ,gi,92 G E such that gi G Afo, i = 1, 2, 
see Corollary 1.3 of IZ9$. 
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Definition 2.6 A monotone map A : E — > E (possiblly multivalued) is said 
to be m-accretive if for some A > 0, A + XI is surjective (as a graph in 
E x E). 

Remark 2.7 So, A is m-accretive, if and only if for all X strictly positive, 
{I + XA)^ 1 is a contraction on E. 

Now, let us consider the case E = The following is taken from [11] . 

Section 1. 

Proposition 2.8 Let [3 : K — > M. be a monotone (possibly multi-valued) map 
such that the corresponding graph is m-accretive. Suppose that (3(0) = 0. 

Let f G E = L 1 (R). 

1. There is a unique u £ L 1 (IR) for which there is w £ L 1 1 oc (IR) such that 

u - Aw = f in V'(R), w(x) £ -/3(u(x)), for a.e. 

(2.15) 

see Proposition 2 of [11]. 

2. It is then possible to define a (multivalued) operator A := An : E — > E 
where D(A) is the set of u G for which there is w G Lj oc (R) 
such that w(x) G |/3(«(z)) /or a.e. x G K and Aw G L 1 ^). For 
u G -D(^4), we set 

Au = {— — w\w as in the definition of D(A)}. 

This is a consequence of the remarks following Theorem 1 in \11^ . 

In particular, if (3 is single-valued, Au = —^Aj3(u). We will adopt this 
notation also if (3 is multi-valued. 

3. The operator A defined in 2. above is m-accretive on E = L 1 (M.), see 
Proposition 2 of \11^ - 

4- We set Jx = (I + Aj4) _1 , which is a single-valued operator. If f G 
L°°(R), then ||Ja/||oo < \\f\\oo, see Proposition 2 (in) of fTIf . In 
particular, for every positive integer n, ||J"/||oo < ll/lloo- 
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Let us summarize some important results of the theory of non-linear semi- 
groups, see for instance [IHl [U EJ [10] or the more recent monograph [29] . 
which we shall use below. Let A : E — > E be a (possibly multivalued) 
m-accretive operator. We consider the equation 

E u'(t) + A(u(t)), <t <T. (2.16) 

A function u : [0, T] — * E, which is absolutely continuous such that for a.e. 
t, u{t, •) E D(A) and fulfills (12. 16H in the following sense, is called strong 
solution. 

There exists r] : [0, T] — ► E, Bochner integrable, such that r)(t) E A(u{t)) for 
a.e. t E [0, Tl and 



u(t) = u + [ r](s)ds, < t < T. 
Jo 



A weaker notion for (|2.16p is the so-called C - solution, see Chapter IV. 8 
of [29]. In order to introduce it, one first defines the notion of e-solution for 

(gup . 

An e-solution is a discretization 

V = {0 = t < h < . . . < t N = T} 
and an E'-valued step function 



uq : t = t 
uj E D(A) : t e]tj-x,tj] 



u e (t) = 

for which tj — t,-_i < e for 1 < j < JV, and 



0E ^ "J- 1 +Au 3 -,1< j<AT. 
tj tj_i 

We remark that, since A is m-accretive, n e is determined by T> and Uq, see 
Proposition 12.81 1. 

Definition 2.9 A C°- solution of [2J6\) is a function u E C([0,T];E) 
such that for every e > 0, there is an e-solution u £ of A2.16\) with 

\\u(t) - u £ (t)\\ <£, < t < T. 
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Proposition 2.10 Let A be an m-accretive (multivalued) operator on a Ba- 
nach space E. We set again J\ := (7 + \A)~ X , A > 0. Suppose uq G D(A). 
Then: 

1. There is a unique C°- solution u:[0,T]^E of l2~Ml 

2. u(t) = lim 

n-*oc J™ u o uniformly in t G [0,T]. 

n 

Proof. 

1) is stated in Corollary IV. 8. 4. of [29] and 2) is contained in Theorem IV 
8.2 of [29]. ■ 

The notion of C°-(or mild) solution needs to be introduced since the dual E* 
of E = 7y 1 (M) is not uniformly convex. If E* were indeed uniformly convex, 
we could have stayed with strong solutions. In fact, according to Theorem 
IV 7.1 of [29j, for a given no G D {A), there would exist a (strong) solution 
u : [0, T] — > E to (|2,16p . which is a simpler notion to deal with. For the 
comfort of the reader we recall the following properties. 

• A strong solution is a C°-solution, by Proposition 8.2 of |29j . 

• Theorem 1.2 of [T7] says the following. Given u$ £ D(A) and given a 
sequence («q) in D(A) converging to no, then, the sequence of the cor- 
responding strong solutions (u n ) converges to the unique C°-solution 
of the same equation. 

3 A porous media equation with singular coeffi- 
cients 

In this section, we will provide first an existence and uniqueness result for 
solutions to the parabolic deterministic equation (II. ip in the sense of dis- 
tributions for multi-valued m-accretive (3. The proof is partly based on the 
theory of non-linear semigroups, see |11| for the case when (3 is continuous. 

However, the most important result of this section, is an existence and 
uniqueness result for a "non-degenerate" linear equation for measures, see 
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(jl.5p . This technical result will be crucial for identifying the law of the 
process appearing in the probabilistic representation (11. 3p . 

We suppose that f3 has the same properties as those given in the introduction. 
However, (3 is allowed to be multi- valued, hence m-accretive, as a graph, in 
the sense of Definition 12.31 Furthermore, generalizing Assumption 11.11 we 
shall assume the following. 

Assumption 3.1 Let [3 : M. — > 2 K be an m-accretive graph with the property 
that there exists c > such that 

w G (3(u) \w\ < c\u\. (3.1) 

Remark 3.2 In particular, (3(0) = and (3 is continuous at zero. We use 
again the representation [3(u] = & 2 (u)u with $ being a non-negative bounded 
multi-valued map <3? : M — > R. 

Remark 3.3 As mentioned before, if (3 : R —* M is monotone (possibly 
discontinuous), it is possible to complete (3 into a monotone graph. For 
instance, if <&(x) = H(x — e c ), then 



(3(x) 





[0,e c 
x 



x < e c 

X — 6 C 

x > e c 



Since the function [3 is monotone, the corresponding graph is monotone. 
Moreover (3 + id is surjective so that, by definition, [3 is m-accretive. 

Proposition 3.4 Let u G L 1 (R) f] L°°(M) Then there is a unique solution 
in the sense of distributions u G (L 1 f"] L°°)([0, T] x M) of 

{ d t u G \d 2 xx (l3{u)), 
[ u(t,x) = u (x), 

that is, there exists a unique couple (u,rj u ) G ((L 1 f] L°°)([0, T] x M)) 2 such 
that 

u(t, x)<f(x)dx = J uo(x)(p(x)dx + - J ds J x)(p"(x)dx, 

Vip G «S(M) and (3.3) 
Vu(t, x ) S 0(u(t, x)) for dt®dx- a.e. (t, x) G [0, T] x E. 
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Furthermore, Olloo — ll u o||oo for every t G [0,T] and £/iere is a unique 

version of u such that u G C([0,T];L x (R)) (C L x ([0,T\ x R)). 

Remark 3.5 1. We remark that, the uniqueness of u determines the 
uniqueness of n £ (3(u) a.e. In fact, for s,t G [0, T], we have 

V 

rj u (r,-)drj = u(t, ■) -u(s,-), a.e. (3.4) 

Since r\ u G ^ 1 ([0, T] xR), i/iis implies that the function r\ u is dt ® cfe- 
a. e. uniquely determined. Furthermore, since (3(0) = and because (3 
is monotone, for dt®dx a.e. (t,x) G [0,T] xIk /lave 

u(i, x) = =>• r/ u (t, cc) = 

and 

u(t,x)rj u (t,x) > 0. 
///? is continuous then we can take rj u (s,x) = (3{u{s,x)). 

3. This result applies in the Heaviside case where <&{x) = H{x — e c ) and 
in the non- degenerate case <£>(x) = H(x — e c ) + e. 

Proof (of Proposition [33]). 

We first recall that by our assumptions, we have (/? + \id)(R) = R for every 
A > 0. 

1. The first step is to prove the existence of a C°-solution of the evo- 
lution problem (|2.16|) in E = L X (R), with A and D(A) as defined in 
Proposition 12.81 2. Suppose D(A) = L X (R). Then, the existence of a 
C°-solution u G C([0, T];L 1 (M)) is a consequence of Proposition 12.81 3. 
and Proposition 12.101 1. In particular, u belongs to L 1 ([0, T] x R). 

2. We now prove that D(A) is dense in in E = L X (R). 

Let u G E. We have to show the existence of a sequence (u n ) in 
D(A) converging to u in E. We set u\ = (I + XA)~ x u, so that u G 
u\ — XA^P(u\). The result follows if we are able to show that 

lim u\ = u, weakly in E, (3.5) 
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because then D(A) is weakly sequentially dense in L 1 (M). In fact we 
can easily show that D(A) is convex and so also its closure. Hence by 
Satz 6.12 of pQ D(A) is also weakly sequentially closed and so the result 
would follow. We continue therefore proving (j3.5[) . Since (/ + A^4) _1 
is a contraction on E, u\ £ E and the sequence (u\) is bounded in 
Since u\ G D(A), by definition, there exists w\ G Lj oc (R) 
such that w\{x) G \[3(u\(x)) for dx-a.e. x G R, Aw\ G L X (R) and 
u = u\ — AAw\. Since (3 has linear growth, w\ also belongs to E for 
every A > and the sequence w\ is bounded in E. Consequently, \w\ 
converges to zero in E when A — > and it follows that \Aw\ converges 
to zero in the sense of distributions, hence u\ — > u again in the sense 
of distributions. Because (u\) is bounded in L 1 (R), it follows that 
u\ — > u weakly in E = L 1 (R), as A — ► 0. 

3. The third step consists in showing that a C°-solution is a solution in 
the sense of distributions of (|3.2|) . 

Let e > and consider a family u £ : [0, T] — > E of e-solutions. Note 
that for «Q := uq and for 1 < j < iV, with ^4 as in Proposition 12.81 2.. 
we recursively have 




(3.6) 



hence 



Awl 




■j-l 



■3-1 



for some w £ j G L 1 1 oc (R) such that Wj G 2/5( n p; dx -a.e. Hence, for 
i G]if_i,if], we have 





o 
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We integrate against a test function a G S(M) and get 



u £ (t,x)a(x)dx = / uo(x)a(x)dx + / / w e (s,x)a"(x)dxds 
Jm. Jo Jm 

(3.7) 

+ (t-tj) / w £ (t £ j ,x)a"(x)dx. 

Letting e go to zero we use the fact that u £ — ► u uniformly in t in 
(u £ ) converges in particular to u G L 1 ([0,T] x R) when e — » 0. 

The third term in the right-hand side of (|3.7j) converges to zero since 



t — t £ is smaller than the mesh e of the subdivision. 



Consequently, (|3.7p implies 

x)a(x)da; = / «o(x)a(x)dx + lim / / tD £ (s, x)a"(x)dxds. 
Jr £ ^°Jo Jm 

(3.8) 

According to our assumption on /3, there is a constant c > such 
that \w £ \ < c|n e |. Therefore the sequence (u; £ ) is equi-integrable on 
[0, T] x R. So, there is a sequence (e n ) such that w £n converges to 
some \r) u G L l ([Q,T} x R) in ^(L 1 , L°°). Taking into account, it 
remains to see that rj u (t,x) G f3(u(t,x)) a.e. ® dx, in order to prove 
that u solves (|3.3p . 

Let K > 0. Using Proposition 12.81 4., by (|3.6p we conclude that 
|| Me (^ ■) II oo < ll^olloo- Consequently for any > 0, the dominated con- 
vergence theorem, implies that the sequence u £n restricted to [0, T] x 
[-K,K] converges to u restricted to [0,T] x [-K, K] in L 2 ([0,T] x 
[—K,K]) and if £n restricted to [0, T] x [— K, K], being bounded by 
c|n e "|, converges (up to a subsequence) weakly in L 2 , necessarily to 
\r] u restricted to [0,T] x [-K,K]. The map v -» i/5(v) on L 2 ([0,T] x 
[— K, K]) is an m-accretive multi- valued map, see [29], p. 164, Exam- 
ple 2c. So it is weakly-strongly closed because of [6j p. 37, Proposition 
1.1 (i) and (ii). Hence the result follows. 

4. The fourth step consists in showing that the obtained solution is in 
L°°([0,T] x R). 
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Point 2. of Proposition [2710] tells us that 



u(t, •) = lim J?uo 

m L 1 (R). Hence, for every t G]0, T] and for some subsequence (n/%) 
depending on i, 

= lim |J^uo| < H^olloo) dx— a.e., 

where we used again Proposition 12,81 4). It follows by Fubini's theorem 
that |u(i,x)| < ||uo||oo> for dt <g> dx-a.e. (t,x) S [0, T] xl. 

5. Finally, uniqueness of the equation in T>'([0, T] x M) follows from The- 
orem 1 and Remark 1.20 of [16J. 



Remark 3.6 1. Theorem 1 and Remark 1.20 of llty apply if (3 is contin- 
uous, to give the uniqueness in point 5. above. However, Remark 1.21 
°f sa D s that this holds true even if 0(0) = and (5 is only contin- 
uous in zero and possibly multi-valued. This case applies for instance 
when <j?(x) = H(x — e c ), e c > 0. 

2. We would like to mention that there are variants of the results in Propo- 
sition \3.4\ known from the literature. However, some of them are just 
for bounded domains while we work in all of BL For instance, when 
the domain is bounded and (3 is continuous, Example 9B in Section 
IV. 9 of i29\j . remarks that a C° -solution is a solution in the sense of 
distributions. 

In order to establish the well-posedness for the related probabilistic rep- 
resentation one needs a uniqueness result for the evolution of probability 
measures. This will be the subject of Theorem 13.81 below. But as will turn 
out, it will require some global L 2 -integrability for the solutions. 

A first step in this direction was Corollary 3.2 of [13], that we quote here for 
the convenience of the reader. 
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Lemma 3.7 Let k g]0, T[. Let p, be a finite Borel measure on [k,T] x R; 
let a, b G L 1 ^, T] x R; p). We suppose that 

I (dtf{t, x) + a(t, x)d xx <p(t, x) + b(t, x)d x <p(t, x)) p(dtdx) = 0, 

J[k,T]xVL 

for all cp G C£°(]0,+oo[xR). Then, there is p G Lf oc ([K,T] x R) such that 

a/ a(t, x)dfj,(t, x) = p(t, x)dtdx. 

We denote the subset of positive measures in .M(R) by .M+(R). 

Theorem 3.8 Let a be a Borel non negative bounded function on [0, T] x R. 
Let Zi : [0,T] — > .M+(R), i = 1,2, 6e continuous with respect to the weak 
topology of finite measures on 

Let z° be an element of .M+(R). Suppose that both z\ and Z2 solve the 
problem dtz = d xx (az) in the sense of distributions with initial condition 
z(0)=z°. 

More precisely, 

/ (p(x)z(t)(dx) = / tf(x)z°(dx) + ds <p" (x)a(s,x)z(s)(dx) (3.9) 
Jr Jr Jo Jr 

for every t G [0,T] and any if G 5(R). 

Then (z± — Z2){t) is identically zero for every t, if z := z\ — Z2, satisfies the 
following: 

ASSUMPTION (A): There is p : [0,T]xR — > R belonging to L 2 ([k, T] xR) 
for every k > such that p(t, •) is the density of z(t) for almost all t g]0, T). 



Remark 3.9 If a > const > 0, then p such that p(t, •) is a density of 
(z\ — Z2){t) for almost all t > 0, always exists, via Lemma \3. 7\ It remains 
to check if it is indeed square integrable on every [k,T] x R. 

Remark 3.10 The weak continuity ofz(t,-) implies that 



SUp ||z(.E)||var < °0- 

te[o,T] 
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Indeed, if this were not true, we could find t n £ [0, T\, such that \\z{t n )\\ Ya3: 
diverges to infinity. We may assume that lirrin-^oo t n = to E [0,T]. Then 

lim / f(x)z(t n )(dx) = / f(x)z(t )(dx) 
n ^°° Jr Jr 

for all f G C&(IR) ; hence by the uniform boundedness principle one gets the 
contradiction that 

SUp || var < OO. 

n 

Remark 3.11 Theorem \3.8\ does not hold without Assumption (A) even in 
the time-homogeneous case. 

To explain this, let $ : M — > R + be continuous and bounded such that <3?(0) = 
and $ is strictly positive on M — {0}. We also suppose that ^ is integrable 
in a neighborhood of zero. 

We choose z° := 5q, i.e. the Dirac measure at zero. It is then possible to 
exhibit two different solutions to the considered problem with initial condition 
z°. 

We justify this in the following lines using a probabilistic representation. Let 
Yq be identically zero. 

According to the Engelbert- Schmidt criterion, see e. g. Theorem 5.4 and 
Remark 5.6 of Chapter 5, ]23^ . it is possible to construct two solutions (in 
law) to the SDE 

Y t = [ $(Y s )dW s . (3.10) 
J o 

where W is a Brownian motion on some filtered probability space. 

One solution is identically zero. The second one Y^ is a non-constant 
martingale starting from zero. We recall the construction of Y^ , since it is 
of independent interest. 

Let B be a classical Brownian motion and we set 




(3.11) 



Problem 6.30 of 123]/ says that the increasing process (Tt) diverges to infinity 
when t goes to infinity. We define pathwise (At) as the inverse of (Tt) and 
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we set Mt = B^ t . M is a martingale since it is a time change of Brownian 

motion. One the one hand we have [M]t = At. But pathwise, by (13.111) we 

have ^ ^ 

A t = I ' <$> 2 (B u )dT u = f <S> 2 (B Av )dv, 
Jo Jo 

through a change of variables u = A v . Consequently we get 

A t = [ <5> 2 {M v )dv. 
Jo 

Theorem 4-2 of Ch. 3 of ]23^ says that there is a Brownian motion W on a 
suitable filtered larger probability space and an adapted process (pt) so that 
M t = J l pdW. We have [M] t = jj p 2 s ds = f* $ 2 (M s )ds, for all t > 0, hence 
p 2 t = $ 2 (Mt) and so $(M t )sign(p t ) = p t - 

We define 

W t = sign(p v )dW v . 
Jo 

Clearly [W]t = t. By Levy's characterization theorem of Brownian motion, 
W is a standard Brownian motion. Moreover, we obtain Mt = Jq Q(M s )dW s 
so that y( 2 ) := M solves the stochastic differential equation kS.luXl . Now 
and Y 2 have not the same marginal laws Vi(t, -),i = 1,2. In fact vi(t, •) is 
equal to 5q for all t £ [0, T] . 

Using Ito's formula it is easy to show that the law v(t, •) of a solution Y of 
113. 1 0\) solves the PDE in Theorem \3. 8\ with a := & 2 and initial condition 5q. 
This constitutes a counterexample to Theorem \3.8\ without Assumption (A). 

Proof (of Theorem 13. 8[ ). 

The arguments developed in this proof is inspired by a uniqueness proof of 
distributional solutions for the porous media equation, see Theorem 1 of |16j . 

Given a locally integrable function (t,x) — > u(t,x), u' (resp. u") stands for 
the first (resp. second) distributional derivative with respect to the second 
variable x. 

In the first part of the proof we do not use Assumption (A). We will explicitly 
state from where it is needed. 
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Let z l ,z 2 be two solutions to (|3.9|) and we set z = z 1 — z 2 . We will study 
the quantity 

g £ (t) = [ B £ z(t)(x)z(t){dx), 



where B £ z(t) G (L 1 P) L°°)(R) is the continuous function v £ defined in Lemma 
T2l taking m = z(t). g e (t) is well-defined, since 



9s(t) < \\z(t)\\ vai sup\B £ z(t)(x)\ for all t G [0,T]. 

Assume we can show that 

lim g £ (t) = for all t G [0, T]. (3.12) 

e^O 

Then we are able to prove that z{t) = for all t G [0, T]. 

Indeed, Lemma 12.21 says that B £ z(t)' is bounded, with a version locally of 
bounded variation and that B £ z(t) G C 6 (K) f| -^ P (M) for all p > 1. 

Let now C,C be positive real constants. Then, since all terms in (I2.11|) are 
signed measures of finite total variation, (|2,lip implies that 



B £ z(t){x)z{t)(dx) = e / {B £ z(t){x)fdx 
c,c] J}-c,c] 



-C,C] 



(3.13) 

B £ z(t)(x)B £ z(t)"(dx). 



If F, G are functions of locally bounded variation, F continuous, G right- 
continuous, classical Lebesgue-Stieltjes calculus implies that 

f FdG = FG(C) - FG{-C) - [ GdF. (3.14) 

J]-C,C] J]-c,c] 

Setting F = B £ z(t),G(x) = B £ z(t)' , we get 

- [ B £ z(t)(x)B £ z(t)"(dx) = -B £ z{t){C)B £ z{t)'{C)+B £ z{t){-C)B £ z{t)'{-C) 
J]-C,C] 

(B £ z{t)'(x)fdx. 



]-C,C] 

Since B £ z(t) G L 1 ^), we can choose sequences (C n ),(C n ) converging to 
infinity such that B £ z(t)(C n ) — ► 0, B £ z(t)(—C n ) — ► as n — > 00. Then, 
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letting n —* oo and using the fact that B £ z(t) and B e z(t)' are bounded, by 
the monotone and Lebesgue dominated convergence theorems, we conclude 
that 

B e z(t)(x)B £ z(t)"(dx) = [ {B £ z{t)'(x)) 2 dx. 



In particular, B £ z{t)' G L 2 (R). Consequently, (|3.13|) implies that 
g&) = j B £ z(t)(x)z(t){dx) =e J(B £ z(t){x)) 2 dx 
+ I {B £ z{t)'{x)fdx. 



In particular, the left-hand side is positive. Therefore, if for all t G [0,T], 
g £ {t) — > 0, as e — ► 0, then 

s/eB £ z(t) -> 

5 e z(i)' -> 

in L 2 (R), as e -> 0, and so, for all t G [0, T], 

z(t) = e£ £ z(i) - B £ z(t)" -» 

in the sense of distributions. Therefore, z = 0. 
It remains to prove (|3.12p . 

Let (5 > and <j>g G (JH) , 0$ > 0, symmetric, with f R cj)$(x)dx = 1 weakly 
approximating the Dirac-measure with mass in x = 0. Set 

z s (t,x):=(<l>s*z(t))(x):= [ <f>s{x - y)z(t)(dy), xeR,te[0,T]. 

JR 

We define h : [0, T] yVf(R) by h(t)(dx) = a(t, x)z(t, dx). Note that by 
(133]) . since 'M^ - •) G <S(R), Vz G R, we have 



z 5 (t,x)= / / $(x-y)h(s)(dy)ds = [ {$*h(s))(x)ds, Vt G [0, T], z G R, 
jo Jr Jo 

(3.15) 

where we used that 21,5(0) = 0, because z(0) = 0, and that x i— > Z5(i,x) is 
continuous for all i G [0, T]. In fact, one can easily prove that z$ is continuous 
and bounded on [0, T] x R. 

Let us now consider u; G <S(R). By Fubini's theorem, for all t G [0, T] it 
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follows that 

w(x)B £ z(t)(x)dx = I w(x) I K £ {x — y)z(t)(dy)dx 



(w*K £ ){y)z{t)(dy). 

Now, B £ z(0) = since z(0) = 0. Therefore by and the fact that 

w * K £ £ S(M), the previous expression is equal to 

/ (w * K £ )" (y)h(s)(dy)ds = I I w" (x)B £ h(s)(x)dxds, 
o Jr Jo Jr 

which in turn by Lemma 12.21 is equal to 

ft r 

w(x)(eB £ h(s)(x)dx — h(s)(dx))ds. 
Consequently, by approximation, 

w(x)B £ z(t)(x)dx = / / w{x){eB £ h{s){x)dx — h(s){dx))ds 

Jo Jr 

(3.16) 

V we C b (R),te [0,T]. 

As a consequence of (I3.15P and (13. 16h and again using Fubini's theorem, for 
all t £ [0, T] we obtain 



o 



g e ,s(t) ■= / z s (t,x)B £ z(t)(x)dx 



arm 







z$(t, x)(eB £ h(s)(x)dx — h(s)(dx))ds 



zs(s,x)(eB £ h(s)(x)dx — h(s)(dx))ds 

Km, 1 " 

+ f I I {<p'^h{r)){x)dr{eB £ h{s){x)dx-h{s){dx))ds 



o 







zs(s,x)(eB £ h(s)(x)dx — h(s)(dx))ds 



pt pr 

+ I {(j>'s*h(r))(x)(eB £ h(s)(x)dx-h(s)(dx))dsdr 

Jo Jo 
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t 

z$(s, x)(eB £ h(s)(x)dx — h{s)(dx))ds 



o 



+ / / {4>s * h(r))(x)B £ z(r)(x)dxdr 



o 

The application of Fubini's theorem above is justified since a is bounded, 
su Pte[o,r] II z (t) II var < °o, K £ is bounded and <f>s G <S(M). But the last term is 
equal to 

t 

(f)'s(x — y)B £ z{r)(x)dxh(r){dy)dr 
4>s(x — y)){eB £ z{r){x)dx — z(r)(dx))h(r)(dy)dr 
£B £ z(r)(x)(4>s * h(r))(x)dxdr — / / z$(r,y)h(r)(dy)dr, 



o 



where we could use Lemma 12.21 in the first step, since 4>s(- — y) £ Vy € 

K. Hence, for all t G [0,T], 

5e,<5(*) = z$(s,x)eB £ h(s)(x)dxds 
Jo Jr 

+ eB £ z(s){x){4>s*h{s))(x)dxds (3.17) 

j JR 

— 2 / z$(s,x)h(s)(dx)ds. 
Jo Jr 

For a signed measure u, we denote its absolute value by \u\. By Lemma 12.21 
we have 

sup / (\z(s)\ * 4>s)(x)eB £ \h(s)\(x)dx < Cy/e, 

s£[0,T] J«. 

where 

C = Trlklloo sup ||z(s)||^ ar , 
z se[o,T] 

and likewise the integrand of the second integral in (|3.17|) is bounded by the 
same constant independent of 5. Hence, as e — * 0, the first and second term 
in the right-hand side of (|3.17p converges to zero uniformly in 5 and uniformly 
in t £ [0, T\. Now, we use Assumption (A), namely that z £ L 2 ([k,T] x R) 
for all k > 0. Then, since B £ z(t) £ L 2 (R),\ft G [k,T], and \\a\loo < oo, (|^TT7|) 
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implies that Vk > 0, t G [k, T] 



9e{t) -ffe(«) 



o— »0 







< 



2^/eTC-2 / I z 2 (s,x)a(s,x)dxds (3.18) 



< 



Now, lim^o 9e( K ) = 0. In fact z(k, •) — ► z(0, •) = weakly, according to the 
assumption of Theorem l3.81 According to Theorem 8.4.10, page 192, of |12j . 
the tensor product z(k, •) <8> z(k, •) converges weakly to zero. On the other 
hand (x,y) i— > K £ (x — y) is bounded and continuous on M 2 . By Fubini's 
theorem 



So, letting first k —>■ in (|3.18p and then e — > 0, (|3.12p follows since ^ e (t) > 
for all t G [0, T]. In fact, we even proved that the convergence in (|3.12p is 



Remark 3.12 Since our coefficient in Theorem 3.6 is only measurable and 
possibly degenerate, to the best of our knowledge this result is really new. 



supposed to satisfy at least Sobolev regularity. 

Theorem 13.81 will be useful for the probabilistic representation of the solution 
of (|3.3p when /3 is non-degenerate. 

4 The probabilistic representation of the deter- 
ministic equation 

Despite the fact that (3 is multi-valued, by its monotonicity and because of 
(|3.1|) . it is still possible to find a multi- valued map $ : R — > R+ such that 

P(u) = $ 2 (u)u, u G R, 

which is bounded, i.e. 




uniformly in t G [0, T]. 



For instance, in recent contributions by f24[ \2G^ . the diffusion coefficient is 



sup sup$(w) < oo. 

«eR* 
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In fact the value of $ at zero is not determined by (5. 

We start with the case where <I> is non-degenerate. The value $(0) being a 
priori arbitrary, we can set 



Definition 4.1 The (possibly) multivalued map (3 (or equivalently is 
called non-degenerate, if there exists some constant cq > such that 
y £ y > Co for any 

4.1 The non-degenerate case 

We suppose in this subsection j3 to be non-degenerate. 

First of all we need to show that solutions of the linear PDE (|3.9p . which 
are laws of solutions to an SDE, are space-time square integrable. 

Proposition 4.2 Suppose a : [0, T] x R — > R to be a bounded measurable 
function which is bounded below on any compact set by a strictly positive 
constant. 

We consider a stochastic process Y = (Y t ,t £ [0, T]) on a stochastic basis 
{Vl,J-,(J-t),P), being a weak solution of the SDE 



Jo 

where W is a standard [Tt)-Brownian motion. For t E [0, T], let z(t) be the 
law ofY t and set z° := z(0). 

1. Then z solves equation \3. 9\) with z° as initial condition. 

2. There is p E L 2 ([0, T] x R) such that p(t,-) is the density of z(t) for 
almost all t E [0, T] . 

3. z is the unique solution of 113. 9\) with initial condition z° having the 
property described in item 2. above. 



$(0) = [liminf inf 3>(u), limsupsup 
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Remark 4.3 A necessary and sufficient condition for the existence and 
uniquess in law of solutions for the equation in Proposition \4-2\ is that Y 
solves the martingale problem of Stroock- Varadhan, see Chap. 6 of ]31^ . re- 
lated to Ltf = a(t,x)f". In our case, existence and uniqueness follow for 
instance from f31]j . Exercises 7.3.2-7.3.4, see a ^ so WB$, Refinements 4-32, 
Chap. 5. We remark that the coefficients are not continuous but only mea- 
surable, so that space dimension 1 is essential. 

The reader can also consult \27\ \28§ for more refined conditions to be able to 
construct a weak solution; however those do not apply in our case. 



Proof (of the Proposition 14. 2\i . 



1. The first point follows from a direct application of Ito's formula to 
tp(Y t ), if G <S(R), cf. the proof of Theorem O 

2. We first suppose that Yq = xo where xo G R. In this case its law 
z° equals 5 X0 i.e. Dirac measure in xq. In Exercise 7.3.3 of |31| . the 
following Krylov type estimate is provided: 



E 



f(t,Y t )dt 



< const 



L 2 ([0,T]xl 



for every smooth function / : [0, T] X R — ► R with compact support. 
This implies the existence of a density (t, y) i— > pt(xo, y) for the measure 
(t,y) » E(f* f(t,Y t )dt). and 

f(t,y)p t (x ,y)dtdy 



L 



[0,T]xl 



< const 



L 2 ([0,T]xM)) 



and const does not depend on xq, but only on lower and upper bounds 
of a. This obviously implies that 



sup 



Pt(x Q ,y)dtdy < oo. 



x eR J [0,T]xR 

This implies assertion 2., when Yq is deterministic. 

If the initial condition Yq is any law z°(dx), then clearly the density of 
Y t is z t (dy) = p(t,y)dy where p(t,y) = f R u {dx)p t (x,y). 
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Consequently, by Jensen's inequality and Pubini's Theorem, 

p 2 (t,y)dtdy < / uo(dx) / pf (t, x,y)dtdy < oo. 

[0,T]xR JR i[0,T]xl 



3. The final assertion follows by 2. from Theorem 137 



Now we come back to the probabilistic representation of equation (II. ip . 

Let us consider the solution u E (L 1 P| L°°)([0, T] x R) from Proposition 13.41 
that is, u solves equation (|3.2p . in the sense of (|3.3p . assuming the initial 
condition uq is an a.e. bounded probability density. Define 



C! if u(£, x) = 0, 

where c\ S ^(0). Note that, because j3 is non-degenerate and Xu(t,x) S 

$(n(t, x)) (8> cix-a.e., we have Xu > Co > 0, (it (8) dx-a.e. Since Xu is only 

defined (it <g> dx-a.e, let us fix a Borel version. According to Remark 14.31 it is 

possible to construct a (unique in law) process Y which is the weak solution 

of t 

Y t = Y + [ Xu(s,Y s )dW s (4.2) 
Jo 

where VP is a classical Brownian motion on some filtered probability space 
and Yo is a random variable so that uq is the density of its law. 

Consider now the law v(t, ■) of the process Y t . We set a(t,x) = Since 
a > c > 0, Proposition 14.21 implies that v £ L 2 ([0,T] x R) and it solves the 
equation 

I dtv = d 2 ,Jav) 
I v(0,x) = u (x). 



On the other hand u itself, which is a solution to (|3.2p (in the sense of (|3.3p ). 
is another solution of equation (14. 3|) . So, being in (L 1 f] L°°)([0, T] x R), u 
is also square integrable. Setting z\ = v,Z2 = u, Theorem 13.81 implies that 
v = u, dt ® dx-a.e. 

Since u G C([0, T], L X (R) and Y has continuous sample paths, it follows that 
u(t, ■) = v(t, •), dx-a.e. for all t E [0, T]. 
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The considerations above prove the existence part of the following represen- 
tation theorem, at least in the non-degenerate case. 



Theorem 4.4 Suppose that Assumption \3. 1\ holds. Let u$ £ L 1 P|L°° such 
that uq > and J R uo(x)dx = 1. Suppose the multi-valued map <J> is bounded 
and non-degenerate. Then there is a process Y , unique in law, such that 
there exists x G (L 1 f] L°°)([0, T] x R) with 

Y t = Y + /„* x (a, y s )^ s (weakly) 

X(i,ar) G *(u(i,ac)), for dt<g>dx-&.e. (t,x) £ [0, T] xl 

Law density of Yt = u(t,-) 

u(0, •) = u , 

(4.4) 

with u G C([0,T];L 1 (IR))ni oo ([0,T] x E). 



Remark 4.5 //$ fs single-valued then \u = 



Proof. Existence has been established above. Concerning uniqueness, 
given two solutions Y\i = 1,2 of (@3D i.e. ([L~3]h By Ui(t,-),i = 1,2, 
we denote the law densities of respectively = 1,2 with corresponding 
Xi and X2- 

The multi- valued version of Theorem II .41 savs that u\ and ti2 solve equation 
(|1.1|) in the sense of distributions, so that by Proposition 13.41 (uniqueness for 
(|3.3p ) we have u\ = U2, and also \i = X2 a.e. 

We note that, since Y^ has a law density for all t > 0, the stochastic integrals 
in (|4,4p are independent of the chosen Borel version of x- Remark 14.31 now 
implies that the laws of Y 1 and Y 2 (on path space) coincide. ■ 

Corollary 4.6 Consider the situation of Theorem \4-4\ and let vo G L 1 f] L°° 
be such that Vq > 0. TTte unique solution v to equation \3. Sty with initial 
condition vq is non negative for any t > 0. Moreover, the mass j K v(t,x)dx 
does not depend on t. 

Proof. Set no = J R vo(y)dy, which we can suppose to be greater than 0. 
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Then the function u(t, x) = solves equation (|3.2p 



Hence, the result follows from Theorem 



d t u = \d 2 xx {^^\ , , 

i 2 xx \ Mo 11 (4 5) 

u(0, ■) = ^ 



Remark 4.7 VFe note i/iai i/ $ is merely bounded below by a strictly posi- 
tive constant on every compact set and if the solutions u are continuous on 
[0, T] x R, then Theorem \4-4\ an d Corollary \4-6\ still hold. In fact, Stroock- 
Varadhan arguments contained in Remark \4-3\ are still valid if Xu is strictly 
positive on each compact set. 



4.2 The degenerate case 

The degenerate case is much more difficult and will be analyzed in detail in 
the forthcoming paper [7J. In this subsection we only explain the first two 
steps in the special case where our j3 of Section 1 is of the form (5{u) = & 2 (u)u 
and the following properties hold: 

Property 4.8 : R — > R is single-valued, continuous on R — {0}. 

Remark 4.9 A priori ^(O) is an interval; however, by convention, in this 
subsection, we will set <3?(0) := liminf £ ^o+ $(u). This implies that <3? is 
always lower semicontinuous. 

We furthermore assume that the initial condition uq and $ are such that we 
have for the corresponding solution u to (jl.ip (in the sense of Proposition 
13. 4|) the following: 

Property 4.10 & 2 (u(t, •)) : R — ► R is Lebesgue almost everywhere continu- 
ous for dt a.e t £ [0, T] . 

Remark 4.11 As will be shown in [1] Property \4- 10[ is fulfilled in many 
interesting cases, for a large class of intial conditions. In fact, we expect to 
be able to show that u(t, •) is even locally of bounded variation if so is uq. 
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Proposition 4.12 Suppose that Property \4-8\ holds. LetuQ > be a bounded 
integrable real function such that f R uo(x)dx = 1 and the corresponding so- 
lution u to |HIP satisfies Property \4-10[ Then, there is at least one process 
Y such that 



Y t = Y + £ $(u(s, Y s ))dW s in law 

Law density(lf) = u(t, ■), 

u(0, •) = u 



(4.6) 



Corollary 4.13 Suppose that Property \4-8\ holds. Letuo £ L 1 f*\L°° be such 
that uq > and that the corresponding solution u to U.l\) satisfies Property 
\4-10 The unique solution u to equation \3. 2\) is non-negative for any t > 0. 
Moreover, the mass f R u(t,x)dx is constant in t £ [0,T]. 



Proof (of Proposition 14. 12|) . We denote the solution to equation (|3.3|) . by 
u = u(t, x). 

Let e e]0,l] and set /3 e (u) = (&(u)+e) 2 u, $ E (u) = &(u)+e,u £ R. Propo- 
sition 13.41 provides the solution u = u £ to the deterministic PDE equation 

f 8 t U = \dl x {(5 £ {u)), 

\ u(0,x) = u (x). 



We consider the unique solution Y = Y £ in law of 

Y t = Y + f^ £ (u(s,Y s ))dW s 

Law density^) = u £ (t,-) 
u 6 (0, •) = u . 



(4.7) 



Since $ + e is non-degenerate, this is possible because of Theorem 14.41 
Since $ is bounded, using Burkholder-Davies-Gundy inequality one obtains 

E{Y t £ - y/} 4 < const(i - s) 2 ,Ve > 0. (4.8) 

where const does not depend on e. Using the Garsia-Rodemich-Rumsey 
lemma (see for instance [8], (3.b), p. 203), we obtain that 



sup E sup 

e>0 \s,ie[0,T] I* 



1*7 -*7I 4 



< oo. 
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Consequently, using Chebyshev's inequality 

lim su P P({ sup \Y t £ - Y s £ \ > A}) = 0, VA > 0. 

5-*0 £> o s,te[0,T},\t-s\<5 

This implies condition (4.7) of Theorem 4.10 in Section 2.4 of [23]. Condition 
(4.6) of the same theorem requires 

lim supP{|y o e | > A} = 0. 

A^+oo £> o 

This is here trivially satisfied since the law of Y £ is the same for all e. Thus 
the same theorem implies that the family of laws of Y £ , e > 0, is tight. 

Consequently, there is a subsequence Y n := Y £n converging in law (as 
C[0, T] -valued random elements) to some process Y. We set <3? n := $ £n 
and u n := u £n where we recall that u n (t,-) is the law of Y™. 

We also set X? = Y t n - Y£. Since 

[X n ] t = f <S>l(u n (s,Y s n ))ds, 
Jo 

and E([X n ]T) is finite, $ being bounded, the continuous local martingales 
X n are indeed martingales. 

By Skorokhod's theorem there is a new probability space (£l,J-,P) and pro- 
cesses Y n , with the same distribution as Y n so that Y n converge P-a.e. to 
some process Y, of course distributed as Y, as C([0,T])- random elements. 
In particular, the processes X n := Y n — Y n remain martingales with respect 
to the nitrations generated by themselves. We denote the sequence Y n (resp. 
Y), again by Y n (resp. Y). 

Remark 4.14 We observe that, for each t G [0, T], u(t, •) is the law density 
ofYf. Indeed, for any t G [0, T], Y™ converges in probability to Y%; on the 
other hand u n (t,-), which is the law of Y™ , converges to u(t,-) in L 1 (M) 
uniformly in t, cf. fT7]/ . Theorem 3 and the preceeding remarks. 

Remark 4.15 Let y n (resp. y) be the canonical filtration associated with 
Y n (resp. Y). 

We set 

f 1 
W n = / - dY n 

1 J $ n ( U n( S ,Y s n)) *■ 
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Those processes W n are standard (y™) - Wiener processes since \W n ]t = t 
and because of Levy's characterization theorem of Brownian motion. Then 
one has 

Yt n = Y£ + I $ n (u n (s,Y s n ))dW?. 
Jo 

We aim to prove first that 

Y t = Y + [ $(u(s,Y 3 ))dW 3 . (4.9) 
Jo 

Once the previous equation is established for the given it, the statement of 
Proposition 14.121 would be completely proven because of Remark 14.141 In 
fact, that Remark shows in particular the third line of (|4.6p . 

We consider the stochastic process X (vanishing at zero) defined by X t = 
Y t — Y . We also set again X? = Y t n - Y n . 

Taking into account Theorem 4.2 of Ch. 3 of [23], as in Remark I3.1H to 
establish (|4.9p it will be enough to prove that X is a y~ martingale with 
quadratic variation [X) t = J & 2 (u(s,Y s ))ds. 

Let s,t € [0, T] with t > s and O a bounded continuous function from 
C([0,s]) to R. 

In order to prove the martingale property for X, we need to show that 

E((X t -X s )G(Y r ,r<s)) = 0. 

But this follows because Y n -> Y a.s. (so X n -> X a.s.) as C([0, T])-valued 
processes; so for each t > 0, X T t l — ► Xt in L 1 (Q) since (X t n , re G N) is bounded 
in L 2 (17) and 

E((X?-X«)@(Y r n ,r<s))=Q. 

It remains to show that X? - f* $ 2 (u(s,Y s ))ds,t £ [0,T], defines a J 7 - 
martingale, that is, we need to verify that 

E ((X? -X 2 S - f $ 2 (n(r, Y r ))dr)Q{Y r , r < s) \ = 0. 
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The left-hand side decomposes into 2(I l (n) + I 2 (n) + I 3 (n)) where 
l\n) = E {{X 2 - X 2 - J <S> 2 (u{r,Y r ))dr)Q(Y r ,r < sfj 

- E (((X?) 2 - (X«) 2 - f\\u(r,Y r n ))dr^ Q(Y r n ,r < s)\ , 
I 2 (n) = I?((V t n ) 2 -ra 2 - jyi(u n (r,Y r n ))dr^Q(Y r n ,r<s)y 



and 



I\n) = E (J* {$l{u n {rX 1 )) ~ <$>\u{rX 1 ))) drQ(Y r n ,r < s) 

We start by showing the convergence of I 3 (n). Now Q(Y™,r < s) is dom- 
inated by a constant. Therefore, since $ n ,^ are uniformly bounded and 
a 2 — b 2 = (a — b)(a + 6), by the Cauchy-Schwarz inequality, it suffices to 
consider the expectation of 

J* ($ n (u n (r, Y r n )) - d>(u(r, Y r n ))) 2 dr (4.10) 

which is equal to 

J* E^ n (u n (r,Y r n )) ~ Hu{r,Y r n ))fdr 

dr I ($ n (u n (r,y)) - $(u(r,y))) 2 u n (r,y)dy. 
Jr 

This equals J\(n) + Jzin) — 2J^(n) where 

Ji(n) = [ dr [ ^ 2 l (u n (r,y))u n {r,y)dy 
Js Jr 

J 2 (n) = I dr I $ 2 (u(r,y))u n (r,y)dy 
Js Jr 

J 3 (n) = dr <S> n (u n (r,y))<S>(u(r,y))u n (r,y)dy. 
Js Jr 



Define 



J := [ [ <S> 2 (u(r,y))u(r,y)dy = [ [ f3{u{r,y))dy 
J s Jr Js Jr 

To show that I 3 (n) — > as n — > oo, it suffices to show that 

lim Ji(n) = lim J 2 (n) = J (4.11) 
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and 

liminf J 3 (n) > J. (4.12) 

Now repeating exactly the same arguments as in Point 3. of Proposition 
E3J it follows that <fr 2 (u n )u n — >■ $ 2 (u)u in ^(L 1 ,^ 00 ) as n — > oo which 
immediately implies (14, lip . 

Furthermore, by Fatou's lemma and since $ n > 

liminf Js{n) > / / liminf 3>(ii n (r, y))$(-u(r, y))u(r, y)dydr 

which by the lower semicontinuity of $>, implies (I4.12|) . 

Now we go on with the analysis of I 2 (n) and I 1 in). I 2 (n) equals zero 
because X n is a martingale with quadratic variation given by [X n ]t = 
j^ 2 n (u n (r,Y r n ))dr. 

We treat finally I^ra). We recall that X n -> X cl. S. cLS cl random ele- 
ment in C([0, T]) and that the sequence E ((X™) 4 ) is bounded, so (X™) 2 are 
uniformly integrable. Therefore, we have 

E ((XD 2 - pC) 2 )9(Y r V < S )) - £ {(X 2 - X 2 )Q(Y r ,r< s)) - 0, 

when n — > oo. It remains to prove that 




E ($ 2 (u(r, y P )) - $ 2 (u(r, Y r n ))e(Y r n , r < s)dr) -» 0. (4.13) 



Now, for fixed dr-a.e. r S [0, T], $(u(r, •)) has a Lebesgue zero set of dis- 
continuities. Moreover, the law of Y r has a density. So, let N(r) be the null 
event of all u E U such that Y r (u) is a point of discontinuity of $(w(r, •)). 
For u ^ iV(r) we have 

lim <5> 2 (u{r,Y r n (uj))) = $ 2 {u(r,Y r (uj))). 

n— »oo 

Hence Lebesgue dominated convergence theorem implies (|4.13p . 



ACKNOWLEDGEMENTS 



36 



The three authors are grateful to the three Referees who gave significant 
suggestions to improve the quality of the paper. 

The second and third named authors would like to thank Prof. Viorel Barbu 
for stimulating discussions. Financial support through the SFB 701 at Biele- 
feld University and nsf-grant 0606615 is gratefully acknowledged. 

References 

[1] H.W. Alt, Lineare Funktionalanalysis: Eine anwendungsorientierte 
Einfiihrung. 4e Auflage. Springer- Ver lag (2002). 

[2] D. G. Aronson, The porous medium equation, in Lect. Notes Math. Vol. 
1224, (A. Fasano and al. editors), Springer, Berlin, 1-46, 1986. 

[3] P. Bak, How Nature Works: The Science of Self- Organized Criticality. 
New York: Copernicus, 1986. 

[4] P. Bantay, I.M. Janosi, Avalanche dynamics from anomalous diffusion. 
Physical review letters 68, no. 13, 2058-2061 (1992). 

[5] V. Barbu, Nonlinear semigroups and differential equations in Banach 
spaces. Noordhoff International Publishing, Leiden, 1976. 

[6] V. Barbu, Analysis and control of nonlinear infinite dimensional systems. 
Academics Press, San Diego, 1993. 

[7] V. Barbu, M. Rockner, F. Russo. Probabilistic representation for solu- 
tions of a porous media equation: the irregular degenerate case. In prepa- 
ration. 

[8] M.T. Barlow, M. Yor, Semimartingale inequalities via the Garsia- 
Rodemich-Rumsey lemma, and applications to local times. J. Funct. Anal. 
49 (1982), no. 2, 198-229. 

[9] S. Benachour, Ph. Chassaing, B. Roynette, P. Vallois, Processu associes 
a V equation des milieux poreux. Ann. Scuola Norm. Sup. Pisa CI. Sci. 
(4) 23, no. 4, 793-832 (1996). 



37 



[10] Ph. Benilan, H. Brezis, M. Crandall, A semilinear equation in L 1 (IR Ar ). 
Ann. Scuola Norm. Sup. Pisa, Serie IV, II Vol. 30, No 2 523-555 (1975). 

[11] Ph. Benilan, M. Crandall, The continuous dependence on (p of solutions 
of u t — Aip(u) = 0. Indiana Univ. Mathematics Journal, Vol. 30, No 2 
161-177 (1981). 

[12] V.I. Bogachev, Measure theory. Vol II. Springer- Verlag 2007. 

[13] V.I. Bogachev, N. Krylov and M. Rockner, On regularity of transition 
probabilities and invariant measures of singular diffusions under minimal 
conditions. Comm. Partial Differential Equations 26 no. 11-12, 2037- 
2080 (2001). 

[14] H. Brezis, Analyse fonctionnelle. Theorie et applications. Masson, 
Paris, 1983. 

[15] H. Brezis, Operateurs maximaux monotones et semigroupes de contrac- 
tion dans les espaces de Hilbert. North Holland, Amsterdam, 1977. 

[16] H. Brezis, M. Crandall, Uniqueness of solutions of the initial-value prob- 
lem for u t - Aip(u) = 0. J. Math. Pures Appl. 58, 153-163 (1979). 

[17] M.G. Crandall, L.C. Evans, On the relation of the operator + to 
evolution governed by accretive operators. Israel Journal of Mathematics 
Vol. 21, No 4, 261-278 (1975). 

[18] L.C. Evans, Nonlinear evolution equations in an arbitrary Banach space. 
Israel Journal of Mathematics Vol. 26, No 1, 1-42 (1977). 

[19] L.C. Evans, Application of nonlinear semigroup theory to certain partial 
differential equations, M. G. Crandall Ed., Academic Press, NY, pp. 163- 
188, 1978. 

[20] A. Figalli, Existence and uniqueness of martingale solutions for SDEs 
with rough or degenerate coefficients. J. Funct. Anal. 254 no. 1, 109-153 
(2008). 

[21] C. Graham, Th. G. Kurtz, S. Meleard, S., Ph. Protter, M. Pulvirenti, 
D. Talay, Probabilistic models for nonlinear partial differential equations. 



38 



Lectures given at the 1st Session and Summer School held in Montecatini 
Terme, May 22-30, 1995. Edited by Talay and L. Tubaro. Lecture Notes 
in Mathematics, 1627, Springer- Verlag. 

[22] B. Jourdain, Probabilistic approximation for a porous medium equation. 
Stochastic Process. Appl. 89, no. 1, 81-99 (2000). 

[23] I. Karatzas, S.E. Shreve, Brownian motion and calculus. Springer- 
Verlag, Second Edition 1991. 

[24] C. Le Bris, C, Lions, P.-L., Existence and uniqueness of solutions to 
Fokker-Planck type equations with irregular coefficients. Comm. Partial 
Differential Equations 33, no. 7-9, 1272-1317 (2008). 

[25] H.P., Jr. McKean, Propagation of chaos for a class of non-linear 
parabolic equations. Stochastic Differential Equations (Lecture Series in 
Differential Equations, Session 7, Catholic Univ., 1967) pp. 41-57. Air 
Force Office Sci. Res., Arlington, Va. 60.75. 

[26] M. Reed, B. Simon, Methods of modern mathematical physics. II. 
Fourier analysis, self-adjointness. Academic Press, New York-London, 
1975. 

[27] A. Rozkosz, L. Slominski, On weak solutions of one- dimensional SDEs 
with time- dependent coefficients. Stochastics and Stochastics Reports, 
Vol. 42, 199-208 (1993). 

[28] T. Senf, On one- dimensional stochastic differential equations without 
drift and time- dependent diffusion coefficients. Stochastics and Stochas- 
tics Reports, Vol. 43, 199-220 (1993). 

[29] R.E. Showalter, Monotone operators in Banach space and nonlinear 
partial differential equations. Providence, RI: American Math. Soc, 1997 

[30] E.M. Stein, Singular integrals and differentiability properties of func- 
tions. Princeton Mathematical Series, No. 30 Princeton University Press, 
1970. 

[31] D.W. Stroock and S.R.S. Varadhan, Multidimensional Diffusion Pro- 
cesses. Springer- Verlag, 1979. 



39 



[32] A.-S. Sznitman, Topics in propagation of chaos. Ecole d' ete de Prob- 
abilities de Saint-Flour XIX— 1989, 165-251, Lecture Notes in Math., 
1464, Springer, Berlin, 1991. 

[33] H. Triebel, Interpolation Theory, Function Spaces, Differential Opera- 
tors. North Holland, Amsterdam 1978. 



40 



